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1 Introduction 


In quantum mechanics, multiple dimensional or multiple particles systems are characterized by the tensor 
product of the Hilbert subspaces [I], where each subspace is associated to each element. It is well known [5] 
that when there is not interaction among these elements, the wave function is just the the tensor product 
of the wave functions of each element, that is, the tensor product of the wave function associated to each 
element determines the non interacting characteristic of the elements in a quantum system. If interaction 
occurs at some time among these elements , this tensor product disappears, and the wave function becomes 
entangled [3]. So, it is necessary to point out that if the wave function is not factorized , the wave function is 
entangled. In this way, the characterization of factorization is somewhat equivalent to the characterization of 
entanglement. In this paper, we will follow this line of ideas to determine the characterization of an entangled 

state [1 El mill m mini HI]- 


2 Factorized State 


When a quantum systems is made up of several quantum subsystems, where the ith-subsystem is characterized 
by a Hamiltonian Hi and a Hilbert space fi, corresponding a two-states and having a basis {|Ci)}{i=o,ii the 
Hilbert space is written as the tensorial product of each subsystem, £ = £n ® ® £i (for n-subsystems). 

The state in each subsystems is defined as a qubit in quantum computation and information theory m and 
is given by 

IV'i) = fli|0) + 6i|l), \0'i\^ + \^i\^ = Ij (1) 

where {|0), |1)} is the basis of the two-states Hilbert subspace £i. A general state Idt) in the Hilbert space £ 
can be written as 

^ith ^|Qp = l, (2) 

where C'^s are complex numbers, and |^) is an element of the basis of £, 

IO = ICn...ei) = len)®---®|ei), efc= 0 ,l fc = l,...,n. ( 3 ) 


A full factorized state in this space is 
where |'0fc) is given by ([1]). 


1 ^') = IV’n) ®---® I'i/'l), 


(4) 


For n = 2, one has a general state jdt) G 

l^r) = CilOO) + C2IOI) + C3IIO) + C4III). (5) 

Let us assume that this state can be written as 

|4') = IV' 2 ) ® IV’i), (6) 

with = 1,2 given by ([T]). Then, substituting ([T]) in ([6]) and equaling coefficients with (0, it follows 

that 

Cl = 0102, C 2 = Oi&2) C^ = & 1 O 2 , Ci = 6i&2- (7) 

From these expression one obtains a single condition for factorization, 

CiCi = C2C3. (8) 

Thus, if this condition is not satisfied, the state (EJ represents an entangled state. So, one can use the 
following known expression [26) as a characterization of an entangled state 

4') =2|CiC4-C2C3|. (9) 
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For n = 3, a general state in the Hilbert space £ is 


l^f) = Ci|00) + C 2 IOOI) + CalOlO) + C'4|011) + CsIlOO) + CellOl) + CrlllO) + Csllll). (10) 


Assuming that this wave function can be written as |5') = ® \ip 2 ) ® IV’i) with ipk) given by ([TJ, and after 

some identifications (as before) and rearrangements, one gets 


C2 Ca Ce Cs h ’ 

C3 C4 C7 Cs 62 ’ 

Cs Cs C7 Cs h ’ 


63 0 

62 0 

61 ^0 


(lla) 

(llb) 

(llc) 


These expression reflex a parallelism between the complex vectors (Ci, Cs, Cs, Cy) and ((72, (74, Cq, Cs), the 
vectors {Ci,C 2 ,C 5 ,Ce) and {C3,C4,C7,Cs), and the vectors ((7i, (72, (73, ( 74 ) and {C 5 ,Cs,C 7 ,Cs)- In addi¬ 
tion, they bring about the following eight independent relations 


(7i(74 - C2C3 = 0 , ( 7 i (76 - (72^5 = 0, Ci (78 - (72^7 = 0, C^Cs - C4C5 = 0 (12) 

(73(78 - (74^7 = 0 , C5Cs-CeCr = 0, Ci(77 - (73C5 = 0 , ^2(78 - (74^6 = 0 . (13) 

If one of these expression is not satisfied, the wave function m represents an entangled state. Therefore, 
one can propose the following expression to characterize an entangled state 

4^) = 2IC1C4 - C2(73| + 2|(7i(76 - C2Cs\ + 2|C'iC8 - ^2(77! + 2\C3Ce - (74(75| 

+2\C3Cs - (74(771 + 2|(75(78 - ^6^7! + 2IC1C7 - C3C5I + 2\C2Cs - C4Cs\. (14) 


For n = 4, a general state in the Hilbert space £ is of the form 

= (7i 10000) -f (7210001) -b (73|010) -b (74|0011) 

-b(75|0100) -b (7610101) -b (7710110) -b (7810111) 

-b(79|100) -b (7io|1001) -b (7ii|110) -b (7 i2 |1011) 

-b(7i3|1100) -b (7 i 4|1101) -b (7 i 5|1110) -b (7i6|llll). (15) 

Assuming this function can be expressed as |'I') = \tp 4 ) (g) lips) ® \tp 2 ) <7) IV'i) with \tpk), k = 1, 2 ,3,4 given by 
o, and after some identifications and rearrangements, one gets 

Cl _ C3 _ C5 _ Cj _ Cq _ Cii _ C13 _ Ci5 

C2 C4 Cq Cs CiQ C12 Ci 4 CiG 


Cl _ C2 _ Cq _ Cq _ Cq _ Cio _ C*i3 _ Cl4 

C3 C4 C7 Cs Cii C12 CiQ CiQ 


Cl _ C2 _ Cs _ C4 _ Cq _ CiQ _ Cii _ C12 

Cq Cq C7 Cs Cis Ci 4 CiQ CiQ 


Ox — C2. _ _C3_ _ C4. _ Cx. _ c^ _ C7_ _ Cs 

Cq Cio Ch C12 C13 C14 C15 Cie ’ 
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expressing similar parallelism we mentioned before. Each row gives us 28 relations, having a total of 112 
possible relations, and from these relations, one can get the following 36 independent conditions 


CiCi - C 2 C 3 = 0 C'4C'i3 - C 7 C 10 = 0 

C,Cs - C^Ce = 0 C 4 C 15 - C 3 C 16 = 0 

CiCn - C 3 C 9 = 0 C 5 C 8 - CeCr = 0 

C 1 C 14 - CgCg = 0 C 5 C 15 - C 7 C 13 = 0 

C 2 CS - CiCe = 0 CeCn - C 5 C 12 = 0 

C 2 C 73 - C 5 C 10 = 0 CeCie - CsCu = 0 

C 2 C 16 - CioCs = 0 C 7 C 12 - CsCn = 0 

C 3 C 15 - C 7 C 11 = 0 C 9 C 14 - CioCii = 0 

C 10 C 16 ~ C 12 C 14 = 0 CiiCie — C 12 C 15 = 0 


CiCe - C 2 C 5 = 0 
CiCio - C 2 C 9 = 0 
C 1 C 12 - C 2 C 11 = 0 
C 1 C 15 - C 5 C 11 = 0 
C 2 C 12 — C 4 C 10 = 0 
C 2 C 74 - CeCio = 0 
C3C8 - C4C7 = 0 
C3C13 - C11C5 = 0 
C10C15 - C11C14 = 


C4C14 — CqCi 2 = 0 (16a) 

C 4 C 16 - C 8 C 12 = 0 (16b) 

C5C14 — CqCi 3 = 0 (16c) 

C 5 C 16 - C7C44 = 0 (16d) 

CeC’is ~ C'leC's = 0 (16e) 

C 7 C 16 - CsCiT, = 0 (16f) 

1^9^12 — C'loC'll = 0 (16g) 

C 9 C 15 - C 11 C 13 = 0 (16h) 

0 C 13 C 16 — C 14 C 15 = 0. (16i) 


Again, if one of these expression fail to happen, (1151) will represents an entangled state. Thus, one can propose 
the following expression to characterize an entangled state made up of 4-qubits basis 


4 ^) = 2|CiC4-C2C3|+2|C4Ci3-C7Cio|+2|CiC'6-C2C5|+2|C4Ci4-C6Ci2| 

+ 2 |CiC 8 - C3C6I + 2IC4C15 - C3C16I + 2 |CiCio - C2C9I + 2IC4C16 - C8C12I 
+2|CiCii - C 3 C 9 I + 2 IC 5 C 8 - C 6 C 7 I + 2|CiCi2 - C 2 C 11 I + 2 IC 5 C 14 - CeCisl 
+2|CiCi4 - C^Cel + 2 IC 5 C 15 - C 7 C 13 I + 2|CiCi5 - C 5 C 11 I + 2 IC 5 C 16 - C 7 C 14 I 
+2IC2C8 - C4C6I + 2 |C 6 Cii - C5C11I + 2IC2C12 - C4C10I + 21 ^ 6^15 - CieCsl 
+2IC2C13 - CgCiol + 2 |C 6 Ci 6 - CsCul + 2IC2C14 - CgCiol + 2IC7C16 - CgCisI 
+ 2 IC 2 C 16 - CioCsl + 2\C7Ci2 - CgCnl + 2 IC 3 C 8 - C 4 C 7 I + 2 IC 9 C 12 - CioCnl 
+2IC3C15 - C7C11I + 2IC9C14 - CioCnl + 2 |C' 3 Ci 3 - CnCsl + 2IC9C15 - C11C13I 
+ 2 |CioCi 6 - C12C14I + 2 |CiiCi 6 - C12C15I + 2 |CioCi 5 - CnC^l + 2 \CnCie - CuCi^l- 

(17) 
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As we can see from these examples, the number of conditions needed to characterize a factorized state (or 
entangled state) grows exponentially with the number of qubits. So, characterization of an entangled state 
for n-qubits in general becomes a very hard work. Now, in terms of the density matrix elements, one could 
have the characterization of entangled states made up of 2,3, and 4 qubits as 

P11P44 + P22P33 — 2i?e(pi2P43) (18) 

= 2^ P 11 P 44 + P22P33 — 2i?e(pi2P43) + 2-\/ piipee + P22P55 — 2i?e(pi2P65) 

+ 2\/P 11 P 88 + P22P77 — 2i?e(pi2P87) + 2-\/ P33P66 + P44P55 ” ‘^ReipsiPSb) 

+2\/ P33PS8 + PiiPn — 2i?e(p34P87) + 2-\/ P55P8S + P66P77 — 2i?e(p56P87) 

+ 2\/piiP77 + P33P55 — 2i?e(pi3P75) + 2-\/ P22P&S + P44P66 ” 2i?e(p24P86)- (19) 

Cp ^ = 2\JpiiPii + P22P33 — 2i?e(pi2P43) + 2-\/P44/313,13 + P77Pl0,10 “ 2i?e(p47Pi3uo) 

+ 2\/ P11P66 + P22P55 — 2i?e(pi2P65) + 2 -\/p 44P14,14 + /566P12,12 — 2i?e(p46Pl4,12) 

+ 2-\/ P11P88 + P33P66 — 2i?e(pi3P86) + 2\/ P44dl5,15 + P33dl6,16 “ 2i?e(p43Pi546) I 
+ 2\/ PllPlO.lO + P22P99 — 2i?e(pi2Pl0,9) + 2-\/ P44PI6.I6 + P88Pl2,12 — 2i?e(p48Pl6,12) 

+ 2\/ PllPlipi + P33P99 — 2i?e(pi3Pii_9) + 2-\/p55P874l'p66P77--'2(Re(p56P^ 

+ 2\/ PllPl2,12 + P22Pll,ll — 2i?e(pi2Pl2,ll) + 2-\/ P55P14,14 + P66Pl3,13 ~ 2iie(/556Pl4,13) 

+ 2\/ PllPl4,14 + P99P66 — 2i?e(pi9Pi4^6) + 2-\/ PssPlS.lS + P77Pl3,13 ~ ‘^R^{P57Pl5,13) 

+ 2\/ PllPl5,15 + P55Pll,ll ~ 2i?e(pi5pi5pi) + 2y^ P55Pl6,16 + P77dl4,14 — 2i?e(p57Pi6,14) 

+2\/ P22P88 + P44P66 ~ 2i?e(p24P86) + 2\/p66Pll,ll + P55Pl2,12 — 2i?e(p6,5Pll,12) 

+ 2\/ P22Pl2,12 + P44dl0,10 ~ 2i?e(p24Pl2,lo) + 2\/ P66dl5,15 + Pl6,16P88 ~ 2iie(p6,16Pl5,8) 

+ 2\/ P22Pl3,13 + P55Pl0,10 ~ 2i?e(p25dl3,lo) + 2-\/ P66dl6,16 + P88dl4,14 — 2i?e(p68Pl6,14) 

+ 2\/ P22Pl4,14 + P66dl0,10 “ 2i?e(p2,6Pl4,lo) + 2-\/P77Pl6,16 + P88dl5,15 ~ 2iie(p78Pl6,15) 

+ 2\/ P22Pl6,16 + PIOPOPSS ~ 2i?e(p2,10Pl6,8) + 2-\/ P77Pl2,12 + PSSPll,!! ~ 2iie(p78Pl2,ll) 

+ 2\/ P33P88 + Pi4P77 ~ “^R^iPSiPS?) + 2\/P99Pl2,12 + Pl0,10Pll,ll “ 2i?e(p9poPl2,ll) 
+2\/P33P15,15 + P77dll,ll ~ 2iie(p37Pi54i) + 2Y^p9gpi4_i4 + Pl0,10Pll,ll ~ 2iie(p9poPl4,ll) 

+ 2\/ P33Pl3,13 + Pll,llP55 ~ 2iie(p3pipi3^5) + 2P99Pl5,15 + Pll,lldl3,13 ~ ‘^R^{P9,llPl513) 

+ 2\/Pl0,10Pl6,16 + Pl2,12dl4,14 — 2i?e(pio,l2Pl614) + 2\/Pll,llPl6,16 + Pl2,12Pl5,15 ~ 2i?e(pii42Pl6,15) 
+2\/Pl0,10Pl5,15 + Pll,llPl4,14 — 2i?e(pio,llPl5,14) + 2-\/Pl3,13Pl6,16 + Pl4,14Pl5,15 ~ 2iie(pi3p4pi6,15)- 

(20) 

For other considerations of entanglement multiqubit entanglement see muni HI HI- Figure 1 below shows 
the values of the expressions (HI and m for 50 entangled states made up of 3-qubits basis and with values 
Cj, with j = 1,..., 6 randomly generated. As we can see, the values obtained with the coefficients Cjs and 
with the density matrix elements pnm are the same. Figure 2 shows for the \W) state, 

IVF) =C2|001) + C3|010) + C5|100), |C2p + |C'3p + |C5p = l, (21) 

the possible values of C'(^(*(4'). As we can see, there are four possible maxima corresponding to the values 
Gi = C-i = = ±l/\/3 and four maxima corresponding to the values Cs = 0, Cg = (73 = ±l/-\/2, related 

to semi-factorized state |0) 0 ((72|01) -I- (73|10)). Figure 3 shows the possible values of (7(^(’(4') for the state 

IGFTZ) = (7i|000) + C8|lll), |Gip + |G8p = l. (22) 
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The maximum value of is gotten for Ci = Cg = as one would expect. 



Figure 1: for arbitrary entangled state. 
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Figure 2: for the state \W) = Ci|001) + C 2 IOIO) + CsIlOO) such that 

|Ci|2+|C2p + |C3|2 = l. 
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Figure 3: for the state \GHZ) = CilOOO) + Cgllll) such that |Cip + \Cs\‘^ = 1. 


For a Hilbert space £ generated by n-qubits, defines a continuous function ^ 5R+ with 

5R+ = [0,+oo), c > 1. Since the coefficients of the wave function defines a compact set on the real space 
due to the relation IC'iP = image of this compact set is a compact set in 5R+ [TB]. Thus, a 

normalization factor is possible to introduce on this function to define any compact set [0,c], which it is not 
important. 







3 Dynamical consideration 


Following Lloyd’s idea m, consider a linear chain of nuclear spin one half, separated by some distance and 
inside a magnetic in a direction 6 ( 2 ) = {0,0, Bo{z)), and making and angle 9 with respect this linear 
chain. Choosing this angle such that cos 9 = 1/'/S, the dipole-dipole interaction is canceled, the Larmore’s 
frequency for each spin is different, oJk = 'yBQ{zk) with 7 the gyromagnetic ratio. The magnetic moment of 
the nucleus flk is related with its spin through the relation pk = qSfe, and the interaction energy between 
the magnetic field and magnetic moments is Hint = — ■ B(zfc) = — Xa, If in addition, one has 

first and second neighbor Ising interaction, the Hamiltonian of the system is just [nnn 


H, 


N 

= -Y^ujkS/- 

k^l 


2J 


N-l 

E ^ksu. 


2 J' 


N-2 

E ^kS^k+2, 


(23) 


where N is the number of nuclear spins in the chain (or qubits), J and J' are the coupling constant of the 
nucleus at first and second neighbor. Using the basis of the register of N-qubits, {|^;„,... ,^i)} with = 0,1, 
one has that = (—l)^'“/i|^fc)/2. Therefore, the Hamiltonian is diagonal on this basis, and its eigenvalues 

are 

N N-l , N-2 

- T ^ (24) 

fe=i fe=i fc=i 

Consider now that the environment is characterized by a Hamiltonian Hg and its interacting with the quantum 
system with Hamiltonian Hg. Thus, the total Hamiltonian would he H = Hg + Hg + Hge, where Hge is the 
part of the Hamiltonian which takes into account the interaction system-environment, and the equation one 
would need to solve, in terms of the density matrix, is [20l [2T| 

tH^ = [H,pt], (25) 

where pt = pt{s,e) is the density matrix which depends on the system and environment coordinates. The 
evolution of the system is unitary, but it is not possible to solve this equation due to a lot of degree of 
freedom. Therefore, under some approximations and tracing over the environment coordinates m , it is 
possible to arrive to a Lindblad type of equation ESI 121] for the reduced density matrix p{s) = treipt), 

= [Hg,p] + (26) 

i=l ^ 


where V) are called Kraus’ operators. This equation is not unitary and Markovian (without memory of the 
dynamical process). This equation can be written in the interaction picture, through the transformation 
p = UpW with U = as 

^h^=^p), (27) 

where C{p) is the Lindblad operator 


^(P) = 


E{w; 


- ^V/V.p- -pv,v/ 


(28) 


with V = UVU^. The explicit form of Lindblad operator is determined by the type of environment to consider 
[25] at zero temperature. So, the operators can be Vi = S~ (for dissipation) for the model independent with 
the environment. In this case, each qubit of the chain acts independently with the environment, and one has 
local decoherence of the system. The Lindblad operator is 

1 ^ 
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where and Sf. are the ascend and descend operators such that = US^U'^ = where Clk is 

defined as flk = Wk + + x(‘^fc+2 + ^k-2)" The solutions of the equations are 






Pieit) 


Pnit) 

Pisit) 


Pll(o) + ^ 22 ( 0 ) + ^ 33 ( 0 ) + P44(0) + P55(0) + P66(0) + P77(0) + P88(0) 

— (^55(0) + P66(0) + P77(0) + P88(0))e“'’'^* — (^33(0) + P44(0) + P77(0) + P88 (0))e“'*'^* 
-(P22(0) + P44(0) + P66(0) + P88(0))e-^^‘ + (p77(0) + P88(0))e-(^i+^^)‘ 

+ (P66(0) + P88(0))e-('^1+'^^)* + (p 44(0) + p88(0))e-(^-+^^)‘ - p88(0)e-(^i+^-+^^)‘; 
7iP58(0)e“*'^i^ 


Pl6(0) 

^ 17 ( 0 ) ' P- -T 

\/ii + (j +j') 


'Jll + U + j') 

72/038 (0)e~"'^^° 

y /72 + 4 j 2 

73P28(0)e"*'^i^ 


1 — e 


(dt+i')-7i)i 


o-'|(72+73)i. 


— 5 ( 27 - 72 )*^ 

1 _ g[*(i+i')-73]t 


-5(71+73)*. 


= ( 
,7'l2 V 


g-|( 7 l+ 72 )t. 
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Pl2{t) 


P23{t) 

P25{t) 

P27{t) 

P28{t) 

P33{t) 


P3b{t) 

P36{t) 

Pssit) 

P44{t) 

P4b{t) 

P4&{t) 

P47it) 
P48{t) 
P55{t) 

P58{t) 

Pm{t) 

P67{t) 

P68{t) 

P77{t) 

P78{t) 

P88{t) 


(P22(0) + P44(0) + P66(0) + P88(0))e-'>'^‘ - (p66(0) + P88(0))e-(T'^+^^)‘ 
-(P44(0) + /988(0))e-(^=+^^)‘ + p88(0)e-('^i+^=+'^^)‘; 


P23(0)+ 


7iP67(0)e 

\/ii + (j - i'Y 


_ g(*0-i')-7i)t 


) 


g-5(72+73)t. 


[P25(0) + P47(0) (1 - 

P27(0)e-5(T'i+'>'=+'^3)i. 

P28(0)e-5(T'i+'>'=+2'^^)*; 

(P33(0) + P44(0) + P77(0) + P88(0))e“'’'^* — (p77(0) + P88(0))e“^^^''"'’'^^‘ 
-(P44(0) + p88(0))e-(^=+^^)‘ + p88(0)e-('^^+'^=+^^)*; 


73P46(0)e 

\/il + (j - i'Y ^ ' 


£-5(71+72)*. 


P36(0)e-5(T'i+'^=+'^=')*; 

P38(0)e-i(T'^+2T'=^+^^)*; 

(P44(0) + p88(0))e-(^^+^^)‘ - p88(0)e-(^i+^^+^^)‘; 


P45(0)e-5(T'i+'^=+'^=’)‘; 

P46(0)e-5(T'i+'^=+2T'3)7 

P47(0)e-5(T'^+2T'^+^^)*; 

P48(0)e-5(T'i+2T'^+2T'=*)‘; 

(P55(0) + P66(0) + P77(0) + P88(0))e“'’'^* — (p77(0) + P88(0))e“^^^''"^^^‘ 

-(P66(0) + p88(0))e-(^^+^^)‘ + p88(0)e-('^^+^=+^^)‘; 


P58(0)e-i(27l+72+73)*. 

(P66(0) +P88(0))e-('^i+^^)‘ - p88(0)e-(^i+^^+'^^)‘; 


P 68 ( 0 )e- 5 ( 27 i+72+273)*. 

(P77(0) +P88(0))e-('^i+^^)‘ - p88(0)e-(^^+^^+^^)‘; 

P78(0)e-5(27i+272+73)*. 

P 88 ( 0 )e-(T' 1 +T'=+T'^)‘. 


where (pij are given by 


4>ia = tan 


-1 J+J 


7i 


023 = tan 


-1 J-J 


7i 


016 = tan 


035 = tan 


1 

72 


1 J-J 


017 = tan 


1 J+J 


73 


73 


In our case, we have three qubits space {|C 3 ‘C 2 Ci)}?i=o,i) and our parameter in units 27tMHz are 

wi = 400; W 2 = 200; 013 = 100 J = 10; J' = 0.4 
71 = 0.05; 72 = 0.05; 73 = 0.05. 


the time is normalized by the same factor of 277+1 Hz, and we include in this study the entangle state 

l'f'i) = ^(1000)+ |111) +1001)+ I110)). 


(30) 
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Figure 4 shows the behavior of the entangled states \ W) and \GHZ) as a function of time when this entangled 
state interact with the environment. Purity behavior, tr(p^), is also shown. The system starts as a pure 
entangled state, it evolves in a mixed state and finishes in the pure ground state (|000)), after sharing energy 
with the environment. The states \GHZ) and behave more robust than the state \ W), G^f^ grows since 
other entangled stated contribute to this function. In contrast, starting with the entangled state \ W), there 

(3) 

are not other entangled states which make contribution to the function G), ^ in the dynamics, and one sees 
an exponential decay. 


o 



Figure 4: Gf''^ and Purity for the entangled state \W), \GHZ), and |4'i). 
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4 Conclusions 


We have studied the full factorization of an state made up of up to 4-qubits basic states. We have seen that 
there is an indication that the number of conditions to characterizes a factorized state grows exponentially 
with the number of qubits. For two, three and four basic qubits, we showed the conditions in order to have a 
factorized state, and if any of one of this conditions fails, one gets instead an entangled state. Therefore, an 
entangled state is also characterized by the complement of each of this conditions, and the resulting expression 
has been denoted by (n=2,3,4). This non negative function expressed in terms of the coefficients of the 
wave function or in terms the density matrix elements represents a measurement of the entanglement of any 
wave function made up of basic n-qubits (n = 2, 3,4). Using this fnnction, we study the decay of entangled 
states \ W), \GHZ), and iTi) due to interaction with environment, and we noticed a great different behavior 
of the function Cp , indicating some type of robustness behavior of the states \GHZ) and |4>i). The main 
reason for this different behavior is that the entangled states \GHZ) and |4>i) contain the ground state |000), 
which is the final state in the dynamics. 
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